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We model the cytoskeleton as a fractal network by identifying each segment with a simple Kelvin- 
Voigt element, with a well defined equilibrium length. The final structure retains the elastic char¬ 
acteristics of a solid or a gel, which may support stress, without relaxing. By considering a very 
simple regular self-similar structure of segments in series and in parallel, in 1, 2 or 3 dimensions, we 
are able to express the viscoelasticity of the network as an effective generalised Kelvin-Voigt model 
with a power law spectrum of retardation times, £ ~ r“. We relate the parameter a with the fractal 
dimension of the gel. In some regimes (0 < a < 1), we recover the weak power law behaviours of 
the elastic and viscous moduli with the angular frequencies, G' ~ G" ~ that occur in a variety 
of soft materials, including living cells. In other regimes, we find different power laws for G' and 
G". 


Microrheology measurements on the cytoskeleton of 
the cell revealed interesting weak power law behaviours 
[D (see [H [3] for recent reviews), which are frequently 
associated to the phenomenological “Soft Glassy Materi¬ 
als” (SGM) model [H |3] . Based on the idea of structural 
disorder and metastability, common to all SGMs, this 
model relates the power law exponent of the elastic and 
viscous moduli, G'{w) ^ G"{w) ~ w°', with 0 < a < 1, 
to a mean-field noise temperature x = a-I-1, with a glass 
transition occurring at cc = 1. 

However, rather than a generic fluidic system above 
the glass transition, the cytoskeleton could be more easily 
associated with a polymer network near the sol-gel tran¬ 
sition. Using the ideas of percolation and self-similarity, 
the power law exponent a has been previously related 
to the fractal dimension dy of a flexible polymer cluster. 
This relationship is not unique and depends on the un¬ 
derlying assumptions of the proposed microscopic models 
(monodispersity vs polydispersity, unscreening vs screen¬ 
ing of excluded volume, etc. - see [5] for a review). In 
particular, if a polydisperse polymeric fractal (prescribed 
by bond percolation theory), following Rouse chain dy¬ 
namics (for flexible polymer chains) with fully screened 
hydrodynamic interactions is considered [7], the expo¬ 
nent a can take values between 0 and 1 for df ranging 
from 2.5 to 1.25, respectively. 

More recently, the weak power law behaviour of the 
cytoskeleton has been associated with the “Glassy Worm 
Like Ghain” (GWLG) model [8] (so-called from its anal¬ 
ogy to SGM model). This model defines an average sep¬ 
aration between the crosslinks along the filaments. If the 
relaxation modes have a wavelength shorter than this 
separation, its relaxation time follows the “Worm Like 
Ghain” (WLG) model for semiflexible polymers, which 
largely compose the cytoskeleton. Otherwise, the relax¬ 
ation spectrum is stretched through an effective Boltz¬ 


mann factor with a characteristic energy e that must be 
overcome to induce a conformational change of the net¬ 
work. This model predicts a high frequency regime with 
a power law exponent a = 3/4 (corresponding to WLG 
model) and a low frequency regime with a second power 
law exponent that depends on the phenomenological pa¬ 
rameter e. Extensions of this model deal with the possi¬ 
bility of transient crosslinking between filaments [5]. 

It is known that the cytoskeleton, or a cellular tissue, 
even in equilibrium, supports a certain amount of stress, 
which is imposed by a substrate or other neighbouring 
cells. The structural cytoskeleton filaments must retain 
a solid character, without a full relaxation. In this article, 
we will consider a very simple model of a solid gel, com¬ 
posed of a regular self-similar network of segments with 
well defined lengths and rigid bonds. By identifying each 
segment with a Kelvin-Voigt viscoelastic element (an hy¬ 
pothesis suggested in [TOl m), we are able to express 
the viscoelasticity of the network as an effective gener¬ 
alised Kelvin-Voigt model with a power law spectrum of 
retardation times, £ ~ r“, where a is related to the net¬ 
work power law distribution of lengths (and eventually to 
its fractal dimension), and to the Kelvin-Voigt particular 
characteristics of each segment. This relation is not di¬ 
rect, since in 2 or 3 dimensions we have a large collection 
of Kelvin-Voigt elements in series and in parallel. When 
0 < a < 1, we recover the weak power law behaviours 
G' ~ G" ~ w“. In other regimes, for a < 0 and a > 1, 
we obtain, first analytically (with some approximations), 
and then numerically, different and interesting power law 
behaviours for G' and G". 

Let us initially consider a simple ID self-similar struc¬ 
ture that is defined by first dividing the system size L by 
a number N^. We will get Nl segments of size L/Nl. 
Now, we pick only Kl < Nl segments and repeat the 
procedure n times (see Fig. top, for = 3 and 
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FIG. 1: Top: simple ID self similar structure of segments. 
Bottom: 2D network of segments generated from the ID self 
similar structure represented at the top. 


Kl = 2). We obtain {Nl — segments of size 

li = L/Nl^ (with 1 < i < n — 1 ), and NlK2~^ segments 
of size In = L/N2- 

We may associate each segment, of size I, with a 
Kelvin-Voigt element, composed of a spring of stiffness 
Gi in parallel with a dashpot of viscosity ry;. The (un¬ 
deformed) spring ensures the length 1. The element be¬ 
haves elastically on long times scales and its dynamics 
comes from the viscous element. Its creep compliance 
is Ji{t) = J/(l — 6 “*/"^') [12], where J; = l/G; and 
T; = Vi/Gi is the element retardation time. Because the 
elements are in series, the total creep compliance is just 
the sum of all the elements’ creep compliances [T2|. Tak¬ 
ing into account the repetition of segment sizes, which 
are a consequence of the self similar construction, we may 
write 

n—1 

J{t) = {Nl - Kl) Y. K^Juit) + NLKl-\JiSt) ( 1 ) 

i=l 


This sum may be approximated by an integral by mul¬ 
tiplying it by di = 1. If di is considered to be small, we 
may write 

/■'" di 

Kimj^di ( 2 ) 


The product \K\^di/dl\ corresponds to the density of seg¬ 
ments of size I (per unit size). If we invert the relation 
I = li = we obtain 


1 L di 11 

* ^ IniVi Jl ^ “iniVil 

and we may write 


( 3 ) 

( 4 ) 


Thus, the total creep compliance becomes 

rh 

J{t) - l-^-^Ji{t)dl (5) 


where Iniin — and Imax — ■ 

Let us now suppose, in very general terms, that Ji ~ 
and T = Ti P, with 7 > 0. In these conditions, we have 


J(t)~ J I ^ 

~ ^ (1 - ^ ( 6 ) 


where Tmin = Tn and Tmax = n. 

The function defined by £ ^ ■j-P-ii/'f corresponds to 
the generalised Kelvin-Voigt retardation spectrum [12j . 
If we make the association 


, ID 
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( 7 ) 


we may write £ ^ t“. The rheological response of the 
structure depends essentially on this parameter a, and 
the minimum and maximum retardation times. 

To calculate the elastic and viscous moduli, we deter¬ 
mine first the Laplace transform of the creep compliance: 

/•OO /•Ti 1 

J(s)=/ J{t)e-‘^*dsr^ ^“”^WTT— 

Jo Jt„ s(l + ST) 


Then, the complex creep compliance, which is given by 
the relation J*{w) = {iw)J{i'w). The elastic and viscous 
compliances, defined through J* = J' — iJ'\ are respec¬ 
tively given by 

(9) 

J"W~ (10) 


Finally, the complex modulus is related to the complex 
compliance through G*J* = 1. The elastic and viscous 
moduli, defined by G* = G' -I- iG", may be determined 
from the relations: 


G' = 


J' 

J '2 -H ’ 


G" = 


J" 

J'2 + J"2 


( 11 ) 


For w I/ti {wt I for any r < Tmax = d), the 
integrals of Eq. and 10 are much simplified, and we 
obtain the scalings J' ~ and J" ^ , with J' J". 

The elastic and viscous moduli scale as 


( 12 ) 

For w ^ l/r„ (wt ^ 1 for any r > Tmin = Tn), the 
integrals of Eq. and are again simplified, yielding 
the scalings J' ^ w~^ and J” ^ w~^, with J’ <C J". We 
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obtain then the same scalings for the elastic and viscous 
moduli: 


G’ 


J"2 


G” 



(13) 


The coefficients of G' and G" are dependent of a but 
their scalings are not. In fact, the exponents of w coincide 
with the scalings of a simple Kelvin-Voigt model. 

For 1/ri ^ w ^ I/tVi, we may in some cases extend 
the limits of the integrals of Eq. and 10 to Tmin 0 
and Tmax —>■ oo, which allow us to obtain the results: 


J'{w) ~ “esc™ (if 0 < a < 2) (14) 

J"{w) ~ ^w““sec™ (if — 1 < a < 1) (15) 


When 0 < a < 1, both integrals are well defined and we 
recover the weak power law behaviours: 


G'{w) - G"{w) ^ = tan ™ (16) 

If a < 0, the integral for the elastic compliance (Eq. 

diverges as Tmin = —)■ 0. The elastic compliance is 
then dominated by the smallest retardation times. In this 
case, we have WTn ^ 1, and J' ^ uP. On the contrary, 
if a > 2, then this integral is dominated by the largest 
retardation times, Tmax = ti. We have wti ^ 1, and 
the power law behaviour J' ^ w~‘^. By the same line 


of reasoning, we may determine from Eq. 10 the power 
law behaviours J" ^ for a < —1 and J" ~ 
for a > 1. Applying these results to each interval of a, 
and using the approximations J' 3> J” for a < 0 and 


J' <C J" for a > 1 (which we may infer from Eq. 14 and 


we 

obtain the 

power 

law behaviours: 


G' 


G" - 


{a < —1) 

(17) 

G' 


G" ~ 


(-1 < a < 0) 

(18) 

G' 


G" - 


(0 < a < 1) 

(19) 

G' 


G" ~ 


(1 < a < 2) 

(20) 

G' 


G" - 


(2 < a) 

(21) 


We note that for a < —1 or a > 2, we have a single power 
law behaviour for G' ^ vP and G" ~ w'^, for all values 
of the angular frequency w. Indeed, for these ranges of 
a, the whole structure is entirely dominated by only one 
Kelvin-Voigt element, corresponding respectively to the 
minimum (a < —1) or the maximum (a > 2) retardation 
times. 

The creep compliance of the ID self similar struc¬ 
ture (Eq. corresponds to an effective discrete gen¬ 
eralised Kelvin-Voigt model. The accurate values for the 
exponents of the power law behaviours, G' ^ and 

G” ~ , for any value of a, are shown in Fig. for 
'Tmax/Tmin = 10®. Representative plots of G'{w) and 
G"{w) are shown in Fig. (see m for calculation de¬ 
tails) . 



FIG. 2: Exponents x and y vs a of the power laws G' ~ 
and G” ~ w®', for a generalised Kelvin-Voigt model with 
a power law spectrum of retardation times ~ r“, for 
Tmax/Train = 10®. The highlighted region refers to the weak 
power law behaviour, in which ®«y«a(0<a<l). 
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FIG. 3: Elastic and viscous moduli, G'{w) (solid line) and 
G"{w) (dashed line) for a generalised Kelvin-Voigt model with 
a power law spectrum of retardation times £ ^ t°' , with a = 
—0.5,0.25,0.75,1.5, Tmin = 10“®, Tmax = 10® (in arbitrary 
units). 


To create a 2D or a 3D network structure of segments, 
we may follow the same self similar construction in the 
other dimensions (see 2D network in Fig. bottom). 
The boundaries between the segments are identified with 
the network crosslinks, which allow us to extend the seg¬ 
ments into the interior part of the system, of size or 

The total creep compliance reflects the deformation of 
the 2D or 3D network along the direction of the applied 
force, or stress, which we take to be the direction of the 
ID structure considered before. But now we have an in¬ 
tricate combination of segments in series and in parallel. 

Let us consider first the 2D network structure depicted 
in Fig. bottom. In the first iteration of our construc¬ 
tion, we have the boundary segments of the square of 
size L divided in = 3 parts. At the boundaries be¬ 
tween the segments, we extend new segments into the 
interior of the system, creating a matrix of x AV ad- 
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jacent squares of size L/Nl- We have thus generated 
Ml = Nl + 1 lines in the direction of the force. These 
lines are in parallel. Each line is composed of Nl seg¬ 
ments of size L/Ni^, in series. This regular structure of 
Ml lines in parallel with Nl segments in series, each of 
one with equal creep compliance, gives the matrix creep 
compliance: 


Ji = 1/Gi constant (this choice was done for the rheo¬ 
logical stiffnesses of the SGM model, or the Rouse model 
[13) : in the latter case, it followed from the equipartition 
theorem), i5 = 0 , and we may obtain the interesting result 

= {6 = 0, j = S+1 = 1) (28) 


= ( 22 ) 

At the n**' iteration, after considering all the different 
elements in series and in parallel, we get a surprisingly 
simple result. In fact, the number of segments in series 
cancels with the number of segments in parallel, yielding: 


j2D(t) 


Nl-Kl 
Nl + 1 


71 — 1 




Nl 

Nl + 1 


J+ii) 


(23) 


In the case of the 3D network structure, due to the 
extra dimension, the number of elements in parallel cor¬ 
responds to the square of the number of the elements in 
series. After counting all the contributions at the n*** 
iteration, we have: 


J3°(t) 


Nl-Kl ^^JiM , Nl Jijt) 

(iVL + 1)2 ^ ^ {Nl + 1)2 ^ ^ 


As for the ID case, these sums may also be approxi¬ 
mated by the integrals: 

J Ji{t)^^dlJ 

(25) 

dr 

(26) 


r^{t) 




4 Kl dl 

J V J 


The weak power law behaviours exponents become 


= -, 
7 


a 


3D 
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(27) 


The ID self-similar structure analysed here is simply 
a line of connected segments. However, the 2D and 
3D self-similar structures correspond to networks with 
non-integer fractal dimensions. We may cover the whole 
structures (except a number of lines) with {K\)^ boxes 
{D = 2,3 for the 2D or the 3D network, respectively) 
of size L/N\^, yielding a fractal dimension df = (if 
df > 1). 

The parameters 6 and 7 depend on the choice of our 
particular model. Several possibilities may be consid¬ 
ered. It is reasonable to assume, as in Stokes’s law, that 
rji is proportional to the viscosity of the solvent and to 
the size of the element, rji ~ 1. Then, the retardation 
time Ti = Jir]i ~ leading to 7 = 5 -I- 1. If we take 


There are other reasonable scaling laws. We may for 
instance invoke the idea of springs in series to justify 
Ji ~ 1. In this case, n ^ P, and = (1 + <^//3)/2. 

In this article, we have presented a new paradigm for 
soft solid or gelled materials. We have shown that a 
self similar network, with a power law distribution of 
segment lengths, may lead to a generalised Kelvin-Voigt 
model with a power law spectrum of retardation times, 
£ ~ t“, where a is related to the fractal dimension of 
the network. We recover the weak power law behaviours 
G' ^ G” ^ w°‘ lor 0 < a < I, often observed in the 
cell cytoskeleton, and other soft materials. For a < 0 
or a > I, we also obtain other interesting power law 
behaviours, which are characteristic of this effective gen¬ 
eralised Kelvin-Voigt model. 

The system here presented compares with the “Soft 
Glassy Material” (SGM) model which, in what regards its 
linear viscoelastic regime, may be associated with a gen¬ 
eralised Maxwell model with a power law spectrum of re¬ 
laxation times H ^ 1 (see [IT] for a detailed compar¬ 
ison). With this model, we also recover the weak power 
law behaviours G' ^ G" ~ for 0 < a < I, but we 
have other, substantially different power law behaviours 
for q; < 0 or a > I 13 Ej. The generalised Maxwell 
model (and the SGM model) reflects a more fluidic sys¬ 
tem, which contrasts with the solid like or gelled gener¬ 
alised Kelvin-Voigt model (obtained from a self-similar 
network) described in this article. 

Soft materials, including the cell cytoskeleton, are usu¬ 
ally very complex materials, in which we probably have 
a mixture of microscopic Kelvin-Voigt elements associ¬ 
ated with Maxwell relaxation structures, both in series 
and in parallel, leading to different scaling behaviours in 
different ranges of angular frequencies. Furthermore, the 
cytoskeleton is an active structure from which we may 
expect novel behaviours, at least for particular ranges of 
characteristic times. We hope, however, that this new 
solid like or gelled paradigm model may bring improved 
understanding of these rheological weak power law be¬ 
haviours, that appear so often in so many complex soft 
materials. 
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